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APPARATUS FOR AND METHOD OF MEASURING CLOCK SKEW 

BACKGROUND OF THE INVENTION 

The present invention relates to a clock skew measurement apparatus and 
5 a clock skew measurement method by which a skew between a plurality of clock 
signals distributed by a clock distribution circuit is measured. 

A clock skew has conventionally been estimated statistically using a time 
interval analyzer or a frequency counter. That is, as shown in Fig. 1 , for 
example, when a reference clock signal CLKg from a clock signal source 11 is 

10 distributed and supplied to registers 12j and 12k, each of the clocks under 

measurement CLKj and CLKk in the respective registers 12j and 12k is inputted 
to a time interval analyzer 13, where a timing difference between a zero-crossing 
point of the clock signal under measurement CLKj and a zero-crossing point of 
the clock signal under measurement CLKk is measured, and a fluctuation of the 

1 5 timing difference is measured by a histogram analysis. An example of clock 
skew measurement using a time interval analyzer 13 is described in, for example, 
"Jitter Analysis Clock Solutions", Wavecrest Corp., 1998. 

However, there is a problem in this clock skew measurement method 
using the time interval analyzer that it takes a long time to acquire a number of 

20 data required for the histogram analysis since, it requires an intermediate 

dead-time until clock skew measurement can be repeated. In addition, in the 
clock skew measurement method using the time interval analyzer, a skew 
between clocks having different frequencies cannot be measured. Therefore, a 
new clock skew measurement method is required for a precise control of a local 

25 clock and a global clock. 
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It is an object of the present invention to provide a measurement 
apparatus of a skew between clocks and its method by which a skew between 
clocks can be measured in a short time compared with a conventional apparatus 
and its method. 

5 It is another object of the present invention to provide a measurement 

apparatus of a skew between clocks and its method by which a skew between 
clocks having different frequencies can be estimated. 



SUMMARY OF THE INVENTION 

10 The clock skew measurement apparatus according the present invention 

comprises timing jitter estimation means for estimating timing jitter sequences of 
a plurality of clock signals under measurement, and clock skew estimation 
means to which the plurality of timing jitter sequences are inputted for 
calculating a timing difference sequence between the plurality of timing jitter 

1 5 sequences to output a clock skew sequence. 

In addition, it is desirable that the clock skew measurement apparatus 
includes second clock skew estimation means to which the clock skew sequences 
are inputted for obtaining a difference between the plurality of clock skew 
sequences. 

20 In addition, it is desirable that the clock skew measurement apparatus 

includes frequency multiplication means to which the timing jitter sequence is 
inputted for multiplying a frequency of the timing jitter sequence to output a 
timing jitter sequence of the frequency-multiplied clock signal under 
measurement. 

25 In addition, it is desirable that the clock skew measurement apparatus 



-3- 

includes deterministic clock skew estimation means for estimating a timing error 
between ideal clock edges of the plurality of clock signals under measurement to 
output a deterministic component of clock skew. 

In addition, it is desirable that the clock skew measurement apparatus 
5 includes clock skew detection means to which the clock skew sequence is 
inputted for obtaining clock skew values between the clock signals under 
measurement from the clock skew sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
the timing jitter estimation means comprises analytic signal transformation 

10 means for transforming a clock signal under measurement into a complex 
analytic signal, instantaneous phase estimation means for obtaining an 
instantaneous phase of the analytic signal, linear phase removing means for 
removing a linear instantaneous phase from the instantaneous phase to obtain an 
instantaneous phase noise, and zero-crossing sampling means to which the 

15 instantaneous phase noise is inputted for sampling only the instantaneous phase 
noise data at timings close to zero-crossing timings of a real part of the analytic 
signal to output a timing jitter sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
the analytic signal transformation means comprises bandpass processing means 

20 to which the clock signal under measurement is supplied for extracting only 
components around a fiindamental frequency from the clock signal under 
measurement to band-limit the clock signal under measurement, and Hilbert 
transformation means for Hilbert-transforming an output signal of the bandpass 
processing means to generate a Hilbert pair of the input signal. 

25 In addition, it is desirable in the clock skew measurement apparatus that 
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the analytic signal transformation means comprises time domain to frequency 
domain transformation means to which the clock signal under measurement is 
supplied for transforming the clock signal under measurement into a both-sided 
spectrum signal in frequency domain, band-limiting means for extracting only 
5 components around a positive fundamental frequency of the clock signal under 
measurement in the both-sided spectrum signal, frequency domain to time 
domain transformation means for inverse-transforming an output of the 
band-limiting means into a signal in time domain. 

In addition, it is desirable in the clock skew measurement apparatus that 

10 the analytic signal transformation means comprises a buffer memory to which 
the clock signal under measurement is supplied for storing therein the clock 
signal under measurement, means for extracting the signal in the sequential order 
from the buffer memory such that the signal being extracted is partially 
overlapped with the signal extracted just before, means for multiplying each 

15 extracted partial signal by a window function, means for transforming each 
partial signal multiplied by the window fiinction into a both-sided spectrum 
signal in frequency domain, band-limiting means for extracting only components 
around a positive fiindamental frequency of the clock signal under measurement 
from the transformed both-sided spectrum signal in frequency domain, means for 

20 inverse-transforming an output of the band-limiting means into a signal in time 
domain, and means for multiplying the transformed signal in time domain by an 
inverse number of the window fimction to obtain a band-limited analytic signal. 

In addition, it is desirable that the clock skew measurement apparatus 
includes AD conversion means to which the clock signal under measurement is 

25 inputted for digitizing an analog signal to convert it to a digital signal. 
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In addition, it is desirable that the clock skew measurement apparatus 
includes waveform clipping means to which the clock signal under measurement 
is inputted for removing amplitude modulation components of the clock signal 
under measurement to extract only phase modulation components of the clock 
5 signal under measurement. 

In addition, it is desirable in the clock skew measurement apparatus that 
the clock skew detection means is peak-to-peak detection means for obtaining a 
difference between the maximum value and the minimum value of the supplied 
clock skew sequence 

10 In addition, it is desirable in the clock skew measurement apparatus that 

the clock skew detection means is RMS detection means for obtaining a 
root-mean-square value (RMS value) of the supplied clock skew sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
the clock skew detection means is histogram estimation means for obtaining a 

1 5 histogram of the supplied clock skew sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
a pass bandwidth of the clock signal under measurement is variable in the 
analytic signal transformation means. 

In addition, it is desirable in the clock skew measurement apparatus that 

20 the timing jitter estimation means further includes low frequency phase noise 
removing means to which the instantaneous phase noise is inputted for removing 
low frequency components of the instantaneous phase noise to output the 
instantaneous phase noise from which the low frequency components have been 
removed to the zero-crossing sampling means. 

25 In addition, the clock skew measurement method according to the present 
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invention includes a step of estimating timing jitter sequences of a plurality of 
clock signals under measurement, and a step of calculating a timing difference 
sequence between the plurality of timing jitter sequences to estimate a clock 
skew sequence. 

5 In addition, it is desirable that the clock skew measurement method 

includes a step of obtaining a difference between the plurality of clock skew 
sequences, to estimate a clock skew sequence. 

In addition, it is desirable that the clock skew measurement method 
includes a step of assigning the timing jitter values of the timing jitter sequence 
10 to provide a timing jitter sequence that corresponds to the frequency-multiphed 
clock signal under measurement. 

In addition, it is desirable that the clock skew measurement method 
includes a step of estimating a timing error between ideal clock edges of the 
plurality of clock signals under measurement to estimate a deterministic 
1 5 component of clock skew. 

In addition, it is desirable that the clock skew measurement method 
includes a step of obtaining clock skew values of the clock signals under 
measurement from the clock skew sequence. 

In addition, it is desirable in the clock skew measurement method that the 
20 step of estimating a timing jitter sequence includes a step of transforming a clock 
signal under measurement into a complex analytic signal, a step of obtaining an 
instantaneous phase of the clock signal under measurement from the analytic 
signal, a step of removing a linear instantaneous phase from the instantaneous 
phase to obtain an instantaneous phase noise, and a step of sampling, using the 
25 instantaneous phase noise as an input, only the instantaneous phase noise data at 
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timings close to zero-crossing timings of a real part of the analytic signal to 
output a timing jitter sequence. 

In addition, it is desirable in the clock skew measurement method that the 
step of transforming the clock signal under measurement into an analytic signal 
5 includes a step of extracting only components around a fundamental frequency 
from the clock signal under measurement to band-limit the clock signal under 
measurement, and a step of Hilbert- transforming an output signal of the 
band-limiting processing means to generate a Hilbert pair of the input signal. 

In addition, it is desirable in the clock skew measurement method that the 

10 step of transforming the clock signal under measurement into an analytic signal 
includes a step of transforming the clock signal under measurement into a 
both-sided spectrum signal in frequency domain, a step of extracting only 
components around a positive fundamental frequency of the clock signal under 
measurement in the both-sided spectrum signal, and a step of 

15 inverse-transforming an output of the band-limiting means into a signal in time 
domain. 

In addition, it is desirable in the clock skew measurement method that the 
step of transforming the clock signal under measurement into an analytic signal 
includes a step of storing the clock signal under measurement in a buffer 

20 memory, a step of extracting the signal in the sequential order from the buffer 
memory such that the signal being extracted is partially overlapped with the 
signal extracted just before, a step of multiplying each extracted partial signal by 
a window function, a step of transforming each partial signal multiplied by the 
window function into a both-sided spectrum signal in frequency domain, a step 

25 of extracting only components around a positive fundamental frequency of the 
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clock signal under measurement from the transformed both-sided spectrum 
signal in frequency domain, a step of inverse-transforming the band-limited 
spectrum signal into a signal in time domain, and a step of multiplying the 
transformed signal in time domain by an inverse number of the window fimction 
5 to obtain a band-limited analytic signal. 

In addition, it is desirable in the clock skew measurement method that the 
step of estimating a deterministic component of clock skew between the clock 
signals under measurement includes a step of using the instantaneous linear 
phases of the plurality of clock signals under measurement as inputs and 

10 obtaining a difference between initial phase angles of the instantaneous linear 
phases to obtain a deterministic component of clock skew. 

In addition, it is desirable in the clock skew measurement method that the 
step of estimating a deterministic component of clock skew between the clock 
signals under measurement includes a step of using the timing jitter sequences of 

15 the plurality of clock signals under measurement as inputs and obtaining a 
correlation between the timing jitter sequences to estimate strongly correlated 
clock edges and to obtain an offset value that relates two clock edge with each 
other. 

In addition, it is desirable in the clock skew measurement method that the 
20 step of estimating a deterministic component of clock skew between the clock 
signals under measurement includes a step of using the plurality of clock signals 
under measurement as inputs for obtaining a mean value of zero-crossing timing 
errors between the clock signals under measurement to obtain a deterministic 
component of clock skew. 
25 In addition, it is desirable that the clock skew measurement method 
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includes a step of performing a waveform clipping of the clock signal under 
measurement in order to remove amplitude modulation components of the clock 
signal under measurement and extract only phase modulation components of the 
clock signal under measurement. 
5 In addition, it is desirable in the clock skew measurement method that in 

the step of obtaining the clock skew value , a difference between the maximum 
value and the minimum value of the clock skew sequence is obtained to calculate 
a peak-to-peak value. 

In addition, it is desirable in the clock skew measurement method that in 
10 the step of obtaining the clock skew value , a root-mean-square value of the clock 
skew sequence is obtained to calculate an RMS value. 

In addition, it is desirable in the clock skew measurement method that in 
the step of obtaining the clock skew value , a histogram data of the clock skew 
sequence is obtained. 

1 5 In addition, it is desirable in the clock skew measurement method that the 

step of estimating the timing jitter further includes a step of using the 
instantaneous phase noise as an input for removing low frequency components 
of the instantaneous phase noise. 

The principle of the present invention will be explained. In this 

20 explanation, a clock signal in a microprocessor unit is used as a clock signal 
imder measurement. 



25 



CLOCK SKEW MEASUREMENT METHOD 

First, a clock skew will be defined. As shown in Fig. 2, a clock skew is 
given by the difference between delay times icd^ and Ted*' of the clock signals 
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CLKj and CLKk to arrive at the respective registers 12j and 12k when the clock 
signals CLKj and CLKk are derived from a reference clock signal CLKg of a 
clock signal source 11 as a reference time point of, for example, a clock 
distribution network. That is, the difference is given by an equation (1). 

5 

T4ew(«T) = x^(nT)-xjd(nT) (1) 

Fig. 3 shows by dotted lines a -an ideal reference clock signal CLKg 
having a fundamental period T, and ideal clock signals CLKj and CLKk. Also 

10 shown in Fig. 3 are a difference Ted" (nT) between a -an actual r ising edge of the 
reference clock CLKg and a -an actual r ising edge of the clock signal CLKj, a 
difference icd"^ (nT) between the rising edge of the reference clock CLKg and a 
rising edge of the clock signal CLKk (n=0, 1, 2, •••), and a difference between the 
Ted* (nT) and Xcd'' (nT), i.e., a clock skew Tskew""^ (nT). 

15 As shown in Fig. 4, it is assumed that rising edges of the clock signals 

CLKg, CLKj and CLKk occur at time tcd^ (nT), ted* (nT) and ted*" (nT) respectively, 
and ideal clock edges (clock edges correspond to jitter-free clock signals) of the 
clock signals CLKg, CLKj and CLKk rise at time (nT)g, (nT)j and (nT)k 
respectively. Then the delay time xj (nT) or Xcd'^ (nT) between the clock 

20 departure of the clock signal CLKj or CLKk and the arrival at the respective 
registers 12j or 12k is expressed by following equation. 
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'^idW = t;,(nT)-tf^(nT) 



= kdm-(nT)j 
A<|)j[n] 



^Skew ^ 



t^(nT)-(nT)gJ+{(nT)j-(nT)g} 



A(l)S[n] 



[sec] 



(2) 



T!^H(nT) = t!;H(nT)-t8 (nT) 



= (iiT) - (nT)k ]- [t^ (nT) - (nT)g ]+ {(nT)k - (nT)g } 



^Skew ^ 



l-A(t)S[n] 



2% 



(1 \ 



[sec] 



(3) 



In this case. 



iiL=('>T)j-(nT)g 



[sec] 



(4) 



and 



'iiL=(°T)k-(nT)g 



[sec] 



(5) 



are a time difference between an ideal clock edge of the clock signal CLKj and 
an ideal clock edge of the reference clock signal CLKg, and a time difference 
between an ideal clock edge of the clock signal CLKk and an ideal clock edge of 
the reference clock signal CLKg, respectively, and those time differences 
correspond to deterministic components of clock skew (deterministic clock skew 
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values) that are determined by their respective paths. In addition, A(|)^[n](Tg/27r) 
(=tcd^ (nT)-(nT)g), A(|>^[n](Tj/27r) (=tcd^ (nT)-(nT)j) and A(|)^[n](Tk/27c)(=tc/ 
(nT)-(nT)k) express timing jitter sequences (each unit is second) of the clock 
signals CLKg, CLKj and CLKk, respectively. When the equations (2) and (3) 
are substituted in the equation (1), a clock skew TskeJ'*^ between the clock signal 
CLKj and the clock signal CLKk is expressed by the following equation. 



^Skew ^ 



2k 



-A(|,8[n] 



I ig 1 



^ '^ikew 



A(l)J[n] 



-A<l)S[n] 



^271^ 



[sec] (6) 



-^Skew + 



A(|,k[n] -^^ -A(|)J[n] 



10 The second term 



A^'^ln] 



^^1-A(t.j[n] 
In ) 



_Vj. 



of the equation (6) corresponds to a random spread (random components) of 
clock skew caused by timing jitter of each clock signal. Then in the clock skew 
measurement method according to the present invention, an random distribution 
15 is obtained by combining quantities each representing how much a clock edge of 
each clock signal is deviated from an ideal clock edge, i.e., timing jitters of the 
clock signals. In this case, it is assumed that the fundamental period of the 
distributed clock signal CLKj is generally equal to the fimdamental period of the 
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distributed clock signal CLKk (Tj=Tk). In addition. 



^^L=CT)k-(nT)j 



[sec] 



(7) 



10 



is a difference between a rising edge of an ideal clock of the clock signal CLKj 
and a rising edge of an ideal clock of the clock signal CLKk, and this xskew"'^ is a 
deterministic component of clock skew determined by the paths of clock 
distribution network. 

A deterministic clock skew value TskeW*'^ can be obtained by, for example, 
obtaining instantaneous phases of corresponding two clock signals under 
measurement CLKj and CLKk, and then by calculating a difference between their 
linear phase components. Fundamental components of the clock signals CLKj 
and CLKk are expressed as cosine waves by following equations (8) and (9), 
respectively. 



15 Xj(t) = AjCos((|)-'(t)) = AjCos 



T- 



t+(|)j-A(|)j(t) 



(8) 



(t) = cos((t) (t)) = Aj^ cos 



t + (|)^-A(|)>^(t) 



(9) 



20 



Here, each of instantaneous phases (|)'(t) and (t)^(t) of the respective Xj(t) 
and Xk(t) is expressed by a sum of a linear instantaneous phase component 27it/TL 
(L=j, k) corresponding to a fundamental period Tl (L=j, k), an initial phase angle 
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(t)o^ (L=j, k), and an instantaneous phase noise component A^\t) (L=j, k). 



(l)J(t) = 



t+<|)J-A(|)J(t) 



[rad] 



(10) 



<^'(t): 



t+(|)S-A(t)'^(t) 



[rad] 



(11) 
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An estimation method of an instantaneous phase of a clock signal will be 
explained later. In the equations (10) and (11), when A(|)(t)=0 is assumed, linear 
instantaneous phases of the respective jitter-free clock signals CLKj and CLKk 
can be estimated by following equations. 



'^iinear(t) = 



T- 



[rad] 



(12) 



15 



20 



•l>linear(t) = 



t + ( 



[rad] 



(13) 



In this case, the ideal rising edge of the respective clock signals CLKj and CLKk 
occur at time (nT)j or (nT)k when each linear instantaneous phase of the left side 
becomes (2n7:-7t/2). And the following relations can be obtained from the 
equations (12) and (13), respectively. 
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(nT)j = 



Inn 



[sec] 



(14) 



(nT)k = 



2n7t- 



:YV 
k2%J 



[sec] 



(15) 



Therefore, the equations (14) and (15) are substituted in the equation (7) to 
obtain the deterministic clock skew value iskeJ'*^ by the following equation. 



4L=("T)k-(nT)j 



2n7c- - - 



V 





^^(271; 


^2nj 





[sec] (16) 



In general, the fundamental periods of the respective clock signals CLKj and 
CLKk are the same (Tj=Tk). That is, a deterministic clock skew value between 
the two clock signals under measurement can be obtained as a difference 
between initial phase angles and in the linear instantaneous phases of the 
respective two clock signals under measurement. 

Here, an initial phase angle ^0 of a clock signal under measurement can 
be obtained by applying a linear line fitting by least squares method to an 
instantaneous phase (|)(k) ( k is a digitizing time point), and by selecting so 
that a formula (17) becomes minimum. 
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f;((|,(k)-(a)ok + io)y 
k=l 



(17) 



In this case, an initial phase angle to be obtained is given by the following 
equation (18). 



2N(2N + l)i;(t)(k)-6i:k(t)(k) 
k=i 



■<l>o = 



N(N - 1) 



(18) 



.$0 = 



2(2N + l)2;(^(k)-6Xk(^(k) 

k=l k=l 

N(N - 1) 



(18) 



A parameter estimation by a linear line fitting is described in, for example, 
10 "Random Data: Analysis and Measurement Procedure" by J.S. Bendat and A.G. 

Piersol, 2nd ed., John Wiley & Sons, Inc., p.362, 1986. 

In addition, an initial phase angle of a clock signal under measurement 

x(t) can be obtained by applying a cosine wave fitting by least squares mthod to 

a clock waveform data x(k) or its fimdamental sine wave component, and by 
15 estimating using a maximum likelihood estimation method so that a formula 

(19) becomes minimum. 



N 

Z 

k=lV 



x(k)-Acos 



((2n 



k + ( 



(19) 
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In this case, an initial phase angle to be obtained is given by the following 
equation (20). 



= -arctan 



2x(k)sm— k 
k=l ^ 



N 



271. 



Zx(k)cos— k 
Vk=l ^ . 



(20) 



A parameter estimation by a maximum likelihood estimator is described in, for 
example, "Fundamentals of Statistical Signal Processing: Estimation Theory" by 
S.M. Kay, Prentice-Hall Inc., pp.167-172, 1993. 

In the above description, it is assumed that corresponding clock edges of 

10 the two clock signals under measurement are not spaced apart from each other 
equal to or more than one period. When the corresponding clock edges of the 
two clock signals under measurement are spaced apart from each other equal to 
or more than one period, a deterministic clock skew value is given by a sum of 
the difference between the initial phase angles and an offset time UoffsetTo of a 

15 clock edge, i.e., by the following equation (21). 



"^Skew =(<l'0"<l>0) TT" + "offset To 



[sec] 



(21) 



20 



A clock signal distributed from the clock signal source has a strong 
cause-and-result relationship with a clock signal of the clock signal source. As 
a result, in general, a phase noise (timing jitter sequence) of a distributed clock 
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signal shows a tendency similar to that of a phase noise (timing jitter sequence) 
of the clock signal source. For this reason, timing jitter sequences A^[n] and 
A(t)^[n] of the respective clock signals CLKj and CLKk distributed from the same 
clock signal source show mutually similar tendencies as shown in, for example, 
5 Fig. 5 A and Fig. 5B. Therefore, an amount of offset time noffset of the 

corresponding clock edges of the respective clock signals CLKj and CLKk can 
be estimated by calculating a correlation function between the timing jitter 
sequences A^[n] and A<l)*^[n], and by searching an offset value at which the 
correlation function shows the maximum value. In summary, the amount of 
10 offset noffset of the clock edges can be obtained from an offset position at which 
the cross-correlation between the instantaneous phase noises show the maximum 
value. 

In addition, a deterministic clock skew value can also be obtained by 
obtaining zero-crossing time points of respective clock signals under 

15 measurement, and by calculating a mean value of time differences between 
corresponding zero-crossings. 

Therefore, according to one mode of the clock skew measurement 
method of the present invention, timing jitters A^[n] and Afyn] of the respective 
two clock signals under measurement Xj(t) and Xk(t) shown in Fig. 6 are first 

20 obtained, and then a deterministic clock skew value xskeJ'^ between the two 
clock signals under measurement Xj(t) and Xk(t) is obtained. Next, a timing 
difference between clock edges is obtained by calculating a difference between 
the timing jitter sequences A^[n] and A(|)^[n], and then the random component of 
the clock skew between the clock signals under measurement Xj(t) and Xk(t) is 

25 obtained. A clock skew Tskew ^[n] between the clock signals under 
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measurement is obtained by obtaining a sum of the random component of the 
clock skew and the deterministic component xskew An obtained clock skew 
Tskew''^[n] is shown in Fig. 7. If necessary, an RMS value and a peak-to-peak 
value of clock skew are obtained from the clock skew sequence Tskew^''^[n] . An 
RMS value TskewjiMs^ '' of clock skew is a standard deviation of the clock skew 
Tskew''''[n], and is obtained by the following equation. 



T4ew,RMS = J [^jsj^sL - Tsl^ew ^ [^^^1 (22) 



10 In this case, N is the number of samples of the measured clock skew data, and 
T'skeJ'^ is a mean value. In addition, a peak-to-peak value Tskew,pp^'^ of clock 
skew is a difference between the maximum value and the minimum value of the 
Tskew^'^Ln], and can be obtained by the following equation. 



nL,p? = max„ (T4t [k] ) - min„ (T [k] ) [sec] (23) 

Fig. 8 shows a histogram of clock skew measured by this clock skew 
measurement method. 
20 According to another aspect of the clock skew measurement method of 

the present invention, a clock skew between clock signals each having a different 
frequency from one another can also be measured. Now, in a clock distribution 
network shown in Fig. 9, a system clock signal CLKg from an extemal system 
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clock source 14 is inputted to a clock signal source 11 comprised of a PLL 
(Phase Locked Loop) circuit, where its frequency is multiplied by M. A clock 
signal CLKg whose frequency is M multiple of that of the CLKg is distributed as 
the clock signals CLKj and CLKk to a network, for example, the registers 12j and 
5 12k. Fig. 10a shows the system clock signal CLKq, Fig. 10b shows an ideal 
clock signal that is created by multiplying the frequency of the CLKg by M, and 
Fig. 10c shows the distributed clock signal CLKj that is created by multiplying 
the frequency of the CLKg by M. A0[1] [rad] of the system clock signal CLKg 
represents a timing fluctuation of its edge from an ideal clock edge. Therefore, 

10 the clock shown in Fig. 10b that is created by multiplying the frequency of the 
system clock signal CLKq shown in Fig. 10a by M is, as shown in Fig. 10b, a 
clock of a case in which the system clock signal CLKg is multiplied by 2 in the 
case of illustrated example. Therefore, the number of clock edges is doubled in 
the case of Fig. 10b. In this case, it is sufficient for a rising edge of a newly 

1 5 increased clock to assign the jitter A0[1 ] of a rising edge of the original system 
clock. In the case where the frequency is multiplied by M, one to one 
relationship of A©|"fn/M]] A0[Ln/Mj] and A^[n] can be obtained when the 
A0[1 ] is assigned (M-1 ) times. In this case, -{^fliLircpresents the maximum 
integer that does not exceed x. When a clock skew between the clock signals 

20 CLKj and CLKg is obtained using the equation (6), an equation (24) is obtained. 
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- "^Skew 
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\2n) 





[sec] (24) 
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A deterministic clock skew value Tskew"'' between the clock signals CLKj and 
CLKg is expressed by a time difference between an ideal clock edge (nMT) j of 
the clock signal CLKj and an ideal clock edge (nMT) g of the system clock 
signal CLKg, and can be obtained by the following equation from the initial 
phase angles of the respective clock signals. 



^Skiw=(^T)j-(nMT)G 



-<Po 



2% 



= <P0 



fMTo^ 



27C 



'0 



[sec] 



(25) 



10 In this case, since the clock signal CLKj is a clock signal that is created by 

multiplying the frequency of the system clock signal CLKg by M, a fundamental 
period Tg of the clock signal CLKg is equal to M times of a fundamental period 
Tj of the clock CLKj (TG=MTj). 

In addition, according to further another aspect of the clock skew 

1 5 measurement method of the present invention, a clock skew between the clock 
signals CLKj and CLKk can be measured, using an apparatus that can 
simultaneously measure two channels, by simultaneously sampling only the 
clock signals CLKj and CLKg at first, and then by simultaneously sampling only 
the clock signals CLKk and CLKg. 

20 That is, only the clock signals CLKj and CLKg are simultaneously 

sampled at first, and then a clock skew between the clock signals CLKj and 
CLKg is obtained by the following equation (26) using the second term in the 
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right side of the first equation in the equation (6). 



^Skew 



[n] = T 



Skew 



+ 



A(t,J[n] -L -A(|)S[n] 



[sec] (26) 



10 



Next, only the clock signals CLKk and CLKg are simultaneously sampled 
at first, and then a clock skew between the clock signals CLKk and CLKg is 
similarly obtained by the following equation (27) 



^Skew ['^l ~ '^Skew + 



\2ti) 



A(|)«[n] 



271 



[sec] (27) 



15 Finally, a difference between the clock skew sequences Tskew^"' and Tskew^'^ 
obtained by the equations (26) and (27), respectively is calculated to estimate a 
clock skew between the clock signals CLKj and CLKk using the following 
equation (28). 
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T4ewW=T|4[n]-T|j^[n] 



\2tij 



-A(l)8[n] 



^Skew 



+ 



A(t,J[n] 



In 



-A<l)S[n] 



^271^ 



^( g,k _ gJ ] 
\:Skew ^Skew / 



A(l,'^[n(^l-A(|>8[n] 



^27t^ 



A(l)J[n] 



^271^ 



[sec] 



-A(l)S[n] 
(28) 



As a result, when clock skews among N clock signals are measured, the 
required number of simultaneous samplings is decreased to (N-1) two- charmel. 
simultaneous measurements from nCi (=N(N-l)/2). In addition, this method 
requires only the minimum number of pins when distributed clock signals on a 
semiconductor chip, for example, are taken out to the outside of the chip. 
Therefore, this method is optimum for an evaluation test of a VLSI. In addition, 
the above procedure can be also applied to clock signals each having different 
frequency from one another. 

The clock skew measurement method according to the present invention 
can be used not only for the estimation of a clock skew between distributed 
clocks in a microprocessor unit as described above, but can also be applied to the 
estimation of a clock skew of other signals. 
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TIMING JITTER ESTIMATION METHOD 

Next, a timing jitter estimation method used in the clock skew 
measurement method according to the present invention will be explained. 

A jitter-free clock signal is a square wave having a fundamental 
5 frequency fo. This signal can be decomposed by Fourier analysis to harmonics 
having frequencies fo, 3fo, 5fo, Since a jitter corresponds to a fluctuation of a 
fundamental frequency of a clock signal under measurement, only signal 
components around the fundamental frequency are handled in the jitter analysis. 

A fundamental component of a jittery clock signal (clock signal under 
10 measurement) can be expressed as a sinusoidal wave, assuming that its 
amplitude is A and its fundamental period is To, by the following equation. 



x(t) = Acos((|)(t)) = Acqs 



27t , . , . . ^ 
— t + <|)o-A(|)(t) 



(29) 



15 Here, (|)(t) is an instantaneous phase of the clock signal under measurement, and 
can be expressed by a sum of a linear instantaneous phase component 27tt/To 
containing the fundamental period To, an initial phase angle (|)o (can be made zero 
in the calculation), and an instantaneous phase noise component A(|)(t). 

When the instantaneous phase noise component A(t)(t) is zero, an interval 

20 between rising zero-crossing points of the clock signal under measurement is a 
constant period To. The non-zero A(|)(t) fluctuates the zero-crossing points of 
the clock signal under measurement. That is, a A(|) (nTo) at a zero-crossing 
point nTo represents a time fluctuation of a zero-crossing point, and is called a 
jitter. Therefore, an instantaneous phase noise A(()(t), at zero-crossing points, a 
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timing jitter of the clock signal under measurement can be obtained, by 
estimating an instantaneous phase ^(t) of the clock signal under measurement, 
and by obtaining a difference between an instantaneous phase and a linear phase 
(corresponds to a phase waveform of a jitter-free ideal clock signal) Int/To+^o. 
5 In the timing jitter estimation method used in the present invention, for 

example, a clock signal under measurement x(t) is transformed into a complex 
analytic signal z(t) at first, and an instantaneous phase (|)(t) of the clock signal 
under measurement x(t) is estimated from the analytic signal z(t). A linear line 
fitting by least squares method is applied to the estimated instantaneous phase 

10 waveform data to obtain a linear instantaneous phase (|)iinear(t) corresponding to 
an instantaneous phase waveform of a jitter-free ideal signal. Then a difference 
between the instantaneous phase and the linear instantaneous phase (|)iinear(t) is 
calculated to obtain an instantaneous phase noise A(t)(t) of the clock signal under 
measurement. The instantaneous phase noise waveform A(t)(t) is sampled at 

1 5 timings (approximated zero-crossing points) closest to each zero-crossing point 
of a real part x(t) of the analytic signal z(t) to estimate an instantaneous phase 
noise at a zero-crossing timing nTo, i.e., a timing jitter A(|)[n](=A(()(nTo)). The 
method of obtaining, in this maimer, an instantaneous phase noise A^(t) to 
estimate a timing jitter A(|)[n] was proposed by the inventors of the present 

20 invention, and is described in, for example, "Extraction of Peak-to-Peak and 
RMS Sinusoidal Jitter Using an Analytic Signal Method" by T.J. Yamaguchi, M. 
Soma, M. Ishida, T. Watanabe, and T. Ohmi, Proceedings of 18th IEEE VLSI 
Test Symposium, pp. 395-402, 2000. 

In this timing jitter estimation method, a timing jitter can be estimated 

25 with high accuracy by removing, using waveform clipping means, amplitude 
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modulation (AM) components of the clock signal under measurement in the state 
that phase modulation components corresponding to jitter of the clock signal 
under measurement are kept in the clock signal under measurement. In 
addition, it is desirable to remove low frequency components of the 
instantaneous phase noise using low frequency component removing means. 

INSTANTANEOUS PHASE ESTIMATION METHOD USING ANALYTIC 
SIGNAL 

An analytic signal z(t) of a clock signal under measurement x(t) is 
defined by a complex signal expressed by the following equation (30). 



In this case, j represents an imaginary unit, and an imaginary part x (t) of the 
complex signal z(t) is a Hilbert transform of a real part x(t). 

On the other hand, Hilbert transform of a time waveform x(t) is defined 
by the following equation (31). 



In this case, x (t) is a convolution of the fiinction x(t) and (Vnf). That is, 
Hilbert transform is equivalent to an output at the time when the x(t) is passed 
through a all pass filter. However, the output x (t) in this case has not been 
changed in terms of its spectrum components, but its phase has been shifted by 



z(t) = x(t)+jx(t) 



(30) 




(31) 
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7C/2. 

Analytic signal and Kilbert transform are described in, for example, 
"Probability, Random Variables, and Stochastic Processes" by A. Papoulis, 2nd 
edition, McGraw-Hill Book Company, 1984. 
5 An instantaneous phase waveform (j)(t) of a clock signal under 

measurement x(t) can be obtained from an analytic signal z(t) using the 
following equation (32). 



10 



<|)(t) = tan 



-1 



.x(t). 



(32) 



Next, an algorithm for estimating an instantaneous phase using Hilbert 
transform will be described. First, Hilbert transform is applied to a clock signal 
under measurement x(t). 



15 



x(t) = Acos 



t + (t)o-A<|)(t) 



(33) 



Then a signal x(t) corresponding to an imaginary part of an analytic signal z(t) 
is obtained as follows. 



20 



x(t) = H[x(t)] = Asin 



t + (|)o-A(t)(t) 



(34) 



Then the clock signal under measurement x(t) is transformed into an anal>tic 
signal z(t) as follows. 
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z(t) = x(t) + jx(t) = Acos 



t + (l)o-A(t)(t) 



+ jAsin 



t + (l)o-A(l)(t) 



(35) 



10 



Here, a band-pass filtering process has been applied to the obtained analytic 
signal. This is because a jitter corresponds to a fluctuation of a fiindamental 
frequency of a clock signal under measurement and hence only signal 
components around the fundamental fi'equency are processed in a jitter analysis. 
Next, a phase function (|)(t) is estimated fi^om the obtained analytic signal z(t) 
using the equation (32). 



(l)(t) = 



|?^t+(t)o-A(t)(t) 



mod27t 



[rad] 



(36) 



Here, the (|)(t) is expressed using principal values of phase in the range of -7i to 
+7C, and has a discontinuity point at the proximity of a point where the phase 
1 5 changes fi^om +7c to Finally, by unwrapping (that is, integer multiples of 271 
are appropriately added to principal values (|)(t)) the discontinuous phase function 
(t)(t), a continuous instantaneous phase function (|)(t), fi"om which discontinuity 
has been removed, can be obtained. 



20 



<l>(t) = 



t + (|)o-A<|)(t) 



[rad] 



(37) 
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A phase unwrapping method is described in "The Cephstrum: A Guide to 
Processing" by Donald G. Childers, David P. Skinner, and Robert C. Kemerait, 
Proceedings of IEEE, vol. 65, pp. 1482-1442, 1977. 

5 TRANSFORMATION TO ANALYTIC SIGNAL USING FAST FOURIER 
TRANSFORM 

A transformation from a clock signal under measurement to an analytic 
signal can be realized by a digital signal processing using a transformation of a 
signal in time domain into a signal in frequency domain like Fast Fourier 

10 Transformatio n (FFT) . 

First, FFT is applied to a digitized clock signal under measurement x(t) 
shown in Fig. 1 1 to obtain a both-sided spectrum (having positive and negative 
frequencies) X(f). Fig. 12A shows the obtained both-sided spectrum X(f). 
Next, as shown in Fig. 12B, only data around the ftindamental frequency in the 

1 5 positive frequency components of the spectrum X(f) are retained, and the 
remaining data are replaced by zeros. In addition, the positive frequency 
components are doubled. These processes in frequency domain correspond to 
limiting bandwidth of the clock signal under measurement and transforming the 
clock signal under measurement into an analytic signal Z(f) in time domain. 

20 Finally, by applying inverse FFT to the band-limited signal Z(f), a band-limited 
analytic signal z(t) can be obtained. 

Transformation to an analytic signal using FFT is described in, for 
example, "Random Data: Analysis and Measurement Procedure" by 
J.S. Bendat and A.G. Piersol, 2nd edition, John Wiley & Sons, Inc., 1986. 

25 In addition, when an estimation of instantaneous phase is the object of 
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the process, the process of doubling the positive frequency components can be 
omitted. 

APPROXIMATED ZERO-CROSSING POINT DETECTION METHOD 
Next, an approximated zero-crossing point detection method will be 
5 described. First of all, the maximum value of a real part x(t) of an analytic 
signal of the inputted clock signal under measurement is defined as 100% level, 
and the minimum value is defined as 0% level to calculate 50% level signal 
value VsQo/o as a zero crossing level. A difference between a sample value and 
the 50% level Vsqo/o and a difference between its adjacent sample value and the 

10 50 % level V500/0 of the x(t), i.e., (x(j-l)-V5oo/o) and ( xQyVsox) are calculated, and 
fiirthermore a product of those difference values (xQ-l)-Wso%) x ( x(j)-V5oo/o) is 
calculated. When the x(t) crosses 50% level, i.e., zero-crossing level, the sign 
of its sample value (x(j-l) -V5o%) or (x(j)-V5oo/o) changes from a negative sign to a 
positive sign or from a positive sign to a negative sign. Therefore, when the 

15 product is negative, it is detected that the x(t) has passed the zero-crossing level, 
and a time point j-1 or j at which a smaller absolute value of the sample value 
(x(j-l) -VsQo/o) or (x(j)-V5oo/o) is detected is obtained as an approximated zero 
crossing point. 

20 WAVEFORM CLIPPING 

Waveform clipping means removes AM components from an input signal 
and retains only PM components in the input signal. A waveform clipping is 
performed by: 1) multiplying an analog or digital input signal by a constant, 2) 
replacing a signal value greater than a predetermined threshold value Thl with 

25 the threshold value Thl , and 3) replacing a signal value less than a 
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predetermined threshold value Th2 with the threshold value Th2. Here, it is 
assumed that the threshold value Thl is greater than the threshold value Th2. 

BiaEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing an example of a clock skew measurement 
using a time interval analyzer; 

Fig. 2 is a diagram typically showing a clock distribution network; 

Fig. 3 is a diagram typically showing a timing of clock skew; 

Fig. 4 is a diagram typically showing a relationship between a timing 
jitter and a clock skew; 

Fig. 5A is a diagram showing an example of a timing jitter A(t>^[n] of a 
clock signal under measurement Xj(t); 

Fig. 5B is a diagram showing an example of a timing jitter A^\n] of a 
clock signal under measurement Xk(t); 

Fig. 6 is a diagram showing an example of a clock signal under 
measurement; 

Fig. 7 is a diagram showing an example of a clock skew between clock 
signals under measurement measured by a clock skew measurement method 
according to the present invention; 

Fig. 8 is a diagram showing an example of a histogram of clock skew 
between clock signals under measurement measured by the clock skew 
measurement method according to the present invention; 

Fig. 9 is a diagram typically showing a clock distribution network having 
different clock domains; 

Fig. 10 is a diagram typically showing a principle of clock skew 
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measurement using a frequency multiplication; 

Fig. 1 1 is a diagram showing an example of a digitized clock signal under 
measurement; 

Fig. 12 A is a diagram showing an example of a both-sided power 
5 spectrum of a clock signal under measurement obtained by FFT; 

Fig. 12B is a diagram showing an example of a band-limited single-sided 
power spectrum; 

Fig. 13 is a block diagram showing a functional configuration of an 
embodiment of a clock skew measurement apparatus according to the present 
10 invention; 

Fig. 14 is a flow-chart showing an embodiment of a clock skew 
measurement method according to the present invention; 

Fig. 15 is a block diagram showing a functional configuration of a 
specific example of a deterministic clock skew estimator 102 shown in Fig. 13; 
15 Fig. 16 is a flow-chart showing a processing example of a deterministic 

clock skew estimating step of the step 202 shown in Fig. 14; 

Fig. 17 is a block diagram showing a functional configuration of another 
embodiment of the clock skew measurement apparatus according to the present 
invention; 

20 Fig. 18 is a flow-chart showing another embodiment of the clock skew 

measurement method according to the present invention; 

Fig. 19 is a block diagram showing an example of a functional 

configuration of a timing jitter estimator used in the clock skew measurement 

apparatus according to the present invention; 
25 Fig. 20 is a flow-chart showing an example of a timing jitter estimation 
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method used in the clock skew measurement method according to the present 
invention; 

Fig. 21 is a block diagram showing another example of a functional 
configuration of an analj^ic signal transformer used in the clock skew 
5 measurement apparatus according to the present invention; 

Fig. 22 is a block diagram showing fiirther another example of a 
functional configuration of an analytic signal transformer used in the clock skew 
measurement apparatus according to the present invention; 

Fig. 23 is a flow-chart showing further another example of the analytic 
10 signal transformation method used in the clock skew measurement method 
according to the present invention; 

Fig. 24 is a block diagram showing a portion of another embodiment of 
the apparatus according to the present invention; 

Fig. 25 is a flow-chart showing a portion of another embodiment of the 
1 5 method according to the present invention; and 

Fig. 26 is a block diagram showing another example of a functional 
configuration of the deterministic clock skew estimator 102; 

DESCRIPTION OF THE EMBODIMENTS 

20 Embodiments of the present invention are described below. 

Fig. 13 shows a functional configuration of an embodiment of a clock 
skew measurement apparatus according to the present invention. This clock 
skew measurement apparatus 100 comprises timing jitter estimators 101a and 
101b for estimating timing jitter sequences A^[n] and Afyn] of clock signals 

25 under measurement Xj(t) and Xk(t), respectively, a deterministic clock skew 
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estimator 102 for estimating a timing error between ideal clock edges of the 
respective clock signals under measurement and for estimating a deterministic 
component Tske\v of clock skew, a clock skew estimator 103 to which the timing 
jitter sequences A^[n] and A(|)^[n] are inputted for calculating a timing difference 
5 sequence between those timing jitter sequences A^[n] and A(t)'^[n] to output a 
clock skew sequence TskeJ'^[n], and a clock skew detector 104 for obtaining a- 
clock skew values between the clock signals under measurement from the clock 
skew sequence. In addition, the clock skew detector 104 comprises a 
peak-to-peak detector 105 for obtaining a difference between the maximum 

10 value and the minimum value of the clock skew sequence Tskew''^[n], an RMS 
detector 106 for calculating an RMS value as a standard deviation of the clock 
skew sequence, and a histogram estimator 107 for obtaining a histogram of the 
clock skew sequence. The timing jitter estimators 101a and 101b estimate 
initial phase angles and of the clock signals under measurement Xj(t) and 

15 Xk(t), respectively in addition to the timing jitter sequences A())^[n] and A(|)^[n], 
and output the estimated initial phase angles to the deterministic clock skew 
estimator 102. A specific configuration of the timing jitter estimators 101a and 
101b will be described later. 

Next, the operation in the case where a clock skew between the clock 

20 signals under measurement Xj(t) and Xk(t) is measured using the clock skew 
measurement apparatus 100 of this embodiment will be described. Fig. 14 
shows a processing procedure of the embodiment of the clock skew 
measurement method according to the present invention. First, in step 201, 
initial phase angles V (|)o^ and timing jitter sequences A<l^[n] and Afyn] of 

25 the respective clock signals under measurement Xj(t) and Xk(t) are estimated by 
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the timing jitter estimators 101a and 101b, respectively. Next, in step 202, a 
difference between the initial phase angles ^o' and (j)o^ of the respective clock 
signals under measurement is calculated by the deterministic clock skew 
estimator 102 to estimate a deterministic component iskew''^ of clock skew 
5 between the clock signals under measurement. Next, in step 203, a clock skew 
sequence Tskew'''*^[n] between the clock signals under measurement Xj(t) and Xk(t) 
is estimated by the clock skew estimator 103 from the timing jitter sequences 
A^[n\ and A(|)'^[n], and the deterministic component iskew'''' of clock skew. 
Finally, in step 204, a clock skew value between the clock signals under 

10 measurement Xj(t) and Xk(t) is obtained by the clock skew detector 104 from the 
estimated clock skew sequence TskeW*'^, and the process ends. 

In the step 202 for estimating a deterministic component of clock skew 
between the clock signals under measurement, the deterministic clock skew 
estimator 102 obtains a deterministic component of clock skew between the 

1 5 clock signals under measurement using the equation (16). In addition, in the 
step 202, the deterministic clock skew estimator 102 may obtain, if necessary, an 
absolute value of the equation (16). In addition, in the step 203 for estimating a 
clock skew sequence between the clock signals under measurement, the clock 
skew estimator 103 obtains a clock skew sequence Tskew''''[n] between the clock 

20 signals under measurement using the equation (6). In the step 204 for obtaining 
a clock skew value between the clock signals under measurement, the 
peak-to-peak detector 105 obtains a peak-to-peak value of clock skew using the 
equation (23), the RMS detector 106 obtains an RMS value of clock skew using 
the equation (22), and the histogram estimator 107 obtains a histogram from the 

25 clock skew sequence. Alternatively, an RMS value and/or the peak-to-peak 
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value may be obtained from only the second term of the equation (6). In 
addition, the step 201 for estimating initial phase angles and timing jitter 
sequences of the clock signals under measurement may be replaced by the 
processing procedure shown in Fig. 20. Furthermore, the step 202 for 
5 estimating a deterministic component of clock skew between the clock signals 
under measurement may be replaced by the processing procedure shown in Fig. 
16. 

The clock skew measurement apparatus shown in Fig. 13 can also be 
modified as an apparatus for estimating only random component of clock skew. 

10 In this case, the deterministic clock skew estimator 102 for obtaining a 

deterministic component of clock skew is omitted. Similarly, the clock skew 
measurement method shown in Fig. 14 can also be modified as a method of 
estimating only random component of clock skew. In this case, the step 202 for 
estimating a deterministic component of clock skew from initial phase angles of 

15 the clock signals under measurement is omitted. 

The deterministic clock skew estimator 102 shown in Fig. 13 estimates a 
deterministic component of clock skew from a difference between initial phase 
angles and of the respective clock signals under measurement. However, 
the deterministic clock skew estimator can also be materialized by the 

20 configuration shown in Fig. 15. That is, this deterministic clock skew estimator 
102 to which initial phase angles and and timing jitter sequences A^[n\ 
and A(j)^[n] of the respective clock signals under measurement Xj(t) and Xk(t) are 
inputted comprises an offset estimator 301 for estimating an offset Uoffset between 
corresponding clock edges of the respective clock signals under measurement 

25 from those timing jitter sequences A^[n] and A(|)^[n], and a deterministic clock 
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skew calculator 302 for calculating a deterministic component iskeW*'^ of clock 
skew between the clock signals under measurement from the initial phase angles 
(t)o^ and and the offset Uoffset between the corresponding clock edges estimated 
by the offset estimator 301. 
5 The operation in the case where a deterministic component of clock skew 

between the clock signals under measurement is estimated using this 
deterministic clock skew estimator 102 will be described. Fig. 16 shows the 
processing procedure. First, in step 401, from the inputted timing jitter 
sequences A(j)^[n] and A(l)^[n] of the clock signals under measurement, an offset 

10 position, at which correlation function between those timing jitter sequences 
shows the largest value, is obtained by the offset estimator 301 to estimate an 
offset Uoffset between the corresponding clock edges. Next, in step 402, a 
deterministic component iskew^'^ of clock skew between the clock signals under 
measurement Xj(t) and Xk(t) is calculated by the deterministic clock skew 

1 5 calculator 302 from the inputted initial phase angles and and the offset 
noffset between the corresponding clock edges, and then the process ends. In the 
step 402 for calculating a deterministic component of clock skew between the 
clock signals under measurement, the deterministic clock skew calculator 302 
obtains a deterministic component of clock skew between the clock signals 

20 under measurement using the equation (21). As indicated by dashed lines in 
Fig. 15, instantaneous phase noises A^{t) and A^\t) to be explained in Fig. 19 
may be inputted to the offset estimator 301 to obtain an offset position at which a 
correlative value between those A^(t) and A^\t) becomes the largest value for 
estimating an offset Uoffset. 

25 Fig. 17 shows a functional configuration of another embodiment of the 
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clock skew measurement apparatus according to the present invention. This 
clock skew measurement apparatus 500 comprises timing jitter estimators 101a, 
101b, 101c, and lOld for respectively estimating timing jitter sequences A(|)'[n], 
A(|)^[n], A(|)'^[n], and A(|)^[n] of respective clock signals under measurement Xj(t), 
5 Xg(t), Xk(t), and Xg(t), deterministic clock skew estimators 102a and 102b for 
respectively estimating a timing error Et^'' between ideal clock edges of the clock 
signals under measurement Xj(t) and Xg(t) and a timing error Et^"^ between ideal 
clock edges of the clock signals under measurement xic(t) and Xg(t) to 
respectively estimate deterministic components xskew^'' and xskew^''^ using those 

10 timing errors Et^'' and Et^'*', clock skew estimators 103 a and 103b to which the 
timing jitter sequences A(t)'[n], A<|)^[n], and A(|)''[n], A(|)^[n] are respectively 
inputted for calculating timing difference sequences between the two inputs at 
each clock skew estimator respectively to output clock skew sequences Tskew^'^[n] 
and Tskew^'''[n] respectively, a clock skew estimator 501 to which those clock 

15 skew sequences Tskew^''[n] and Tskew^'''[n] are inputted for obtaining a difference 
between those clock skew sequences to estimate a clock skew sequence 
Tskew''''[n], and a clock skew detector 104 for obtaining a clock skew value 
between the clock signals under measurement from the clock skew sequence 
obtained by the clock skew estimator 501 . For simplicity, the explanation of 

20 portions duplicated with those in Fig. 13 is omitted. 

Next, the operation in the case where a clock skew between the clock 
signals under measurement is measured using the clock skew measurement 
apparatus 500 according to the present invention will be described. Fig. 18 
shows a processing procedure of the clock skew measurement method according 

25 to the present invention. First, in step 601, initial phase angles and ^o^, and 
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timing jitter sequences A^[n] and A(|)^[n] of the respective clock signals under 
measurement Xj(t) and Xg(t) are estimated by the timing jitter estimators 101a and 
101b, respectively. Next, in step 602, a difference between the initial phase 
angles (j)o^ and of the clock signals under measurement is calculated by the 
5 deterministic clock skew estimator 102a to estimate a deterministic component 
'Cskew^"' of clock skew between the clock signals under measurement. Next, in 
step 603, a clock skew sequence Tskew^''[n] between the clock signals under 
measurement is estimated by the clock skew estimator 103a from the timing 
jitter sequences A^[n] and A(t)^[n], and the deterministic component Tskew of 

10 clock skew. 

Next, in step 604, initial phase angles and ^o^, and timing jitter 
sequences A^\r\] and A<|)^[n] of the respective clock signals under measurement 
Xk(t) and Xg(t) are estimated by the timing jitter estimators 101c and 101 d, 
respectively. Next, in step 605, a difference between the initial phase angles 

15 and of the clock signal under measurement obtained in the step 604 is 
calculated by the deterministic clock skew estimator 102b to estimate a 
deterministic component xskew of clock skew between the clock signals under 
measurement. Next, in step 606, a clock skew sequence Tskew^'''[n] between the 
clock signals under measurement is estimated by the clock skew estimator 103b 

20 from the timing jitter sequences A^\n] and A(|)^[n] obtained in the step 604 and 
the deterministic component Tskew ^'^ of clock skew obtained in the step 605. 
Next, in step 607, a clock skew sequence Tskew'^[n] between the clock signals 
under measurement Xj(t) and Xk(t) is estimated by the clock skew estimator 501 
from the clock skew sequences Tskew^''[n] and Tskew^'^'M respectively obtained in 

25 the steps 603 and 606! Finally, in step 608, a clock skew value between the 
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clock signals under measurement Xj(t) and Xk(t) is obtained by the clock skew 
detector 104 from the clock skew sequence TskeJ^^Ln] estimated in the step 607, 
and the process ends. In the step 607 for estimating the clock skew sequence 
between the clock signals under measurement Xj(t) and Xk(t), the clock skew 
5 estimator 501 obtains a clock skew sequence between the clock signals under 
measurement using the equation (28). The process sequence of the steps 
601-603 and the steps 604-606 may be exchanged. For simplicity, the 
explanation of portions duplicated with those in Fig. 14 is omitted. 

The clock skew measurement apparatus shown in Fig. 17 may also be 

10 configured as an apparatus for estimating only random component of clock skew. 
In this case, the deterministic clock skew estimators 102a and 102b for obtaining 
the deterministic components of clock skew are omitted. Similarly, the clock 
skew measurement method shown in Fig. 18 may also be modified as a method 
of estimating only random component of clock skew. In this case, the steps 602 

15 and 605 for estimating deterministic components of clock skew from the initial 
phase angles of the clock signals under measurement are omitted. 

As indicated by dashed lines in Fig. 13, a timing jitter sequence estimated 
by one of the timing jitter estimators 101a and 101b, i.e., the timing jitter 
sequence Afyn] estimated by the timing jitter estimator 101b in the case of Fig. 

20 13 may be copied (M-1) times by the frequency multiplier 701 to obtain a timing 
jitter sequence that is to be obtained when the frequency of the clock signal 
under measurement Xk(t) is multiplied by M, and to supply the obtained timing 
jitter sequence to the clock skew estimator 103. By such a process, in the clock 
distribution system previously described with reference to Fig. 9, the clock 

25 signal Xk(t) in Fig. 13 corresponds to the system clock signal CLKq in Fig. 9, and 
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the clock signal Xj(t) is a signal created by multiplying the frequency of the clock 
signal CLKG(xk(t)) by M. In this case, a clock skew Tskew*^"'[n] between the 
clock signal CLKG(xk(t)) and the clock signal Xj(t) can be obtained by the 
equation (24). 

5 Regarding the processing procedure of the clock skew measurement 

method in this case, as shown in Fig. 14, after timing jitter sequences are 
obtained in the step 201, the timing jitter sequence estimated by the timing jitter 
estimator 101b may be assigned, for example (M-l) times as indicated by dashed 
lines, in step 801 by the frequency multiplier to obtain a timing jitter sequence 

10 that is to be obtained when the frequency of the clock signal under measurement 
is multiplied by M, and then the process may move to the step 202. At this 
time, in the step 202 for estimating a deterministic component of clock skew 
between the clock signals under measurement, the deterministic clock skew 
estimator 102 obtains a deterministic component of clock skew between the 

1 5 clock signals under measurement using the equation (25). In addition, in the 
step 203 for estimating a clock skew sequence between the clock signals under 
measurement, the clock skew estimator 103 obtains a clock skew sequence 
between the clock signals under measurement using the equation (24). 

Also in the case of using the frequency multipher 701, when the clock 

20 skew measurement apparatus is constructed as an apparatus for estimating only 
random component of clock skew, the deterministic clock skew estimator 102 for 
obtaining a deterministic component of clock skew may be omitted. In this 
case, in the clock skew measurement method, the step 202 for estimating a 
deterministic component of clock skew from the initial phase angles of the clock 

25 signals under measurement may be omitted. 
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In addition, the frequency multiplier 701 may also be built in the clock 
skew measurement apparatus shown in Fig. 17. In this case, the frequency 
multiplier is inserted in series in each output side of the timing jitter estimators 
101b and 10 Id. Similarly, a step of multiplying a frequency can also be added 
5 to the clock skew measurement method shown in Fig. 18. In this case, the step 
of multiplying a frequency is inserted after each of the steps 601 and 604 for 
estimating timing jitters. 

Fig. 19 shows an example of configuration of each of the timing jitter 
estimaitors 101a, 101b, 101c, and lOld. This timing jitter estimator 900 is 

10 described in, for example, "Extraction of Peak-to-Peak and RMS Sinusoidal 
Jitter Using an Analytic Signal Method" by T.J. Yamaguchi, M. Soma, M. 
Ishida, T. Watanabe, and T. Ohmi, Proceedings of 18th IEEE VLSI Test 
Symposium, pp. 395-402, 2000. This timing jitter estimator 900 comprises an 
anal5^ic signal transformer 901 for transforming a clock signal under 

15 measurement into a band-limited complex analytic signal, an instantaneous 
phase estimator 902 for obtaining an instantaneous phase of the analytic signal 
transformed by the analytic signal transformer 901, a linear phase remover 903 
for removing a linear instantaneous phase from the instantaneous phase 
estimated by the instantaneous phase estimator 902 to obtain an instantaneous 

20 phase noise, a zero-crossing detector 904 to which a real part of the analytic 
signal is inputted from the analytic signal transformer 901 for generating 
sampling pulses at timings (approximated zero-crossing points) close to 
zero-crossing timings of the real part of the analytic signal, a zero-crossing 
sampler 905 to which the instantaneous phase noise estimated by the linear phase 

25 remover 903 is inputted for sampling the instantaneous phase noise using the 
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sampling pulses from the zero-crossing detector 905 . The analytic signal 
transformer 901 may be constructed such that a pass bandwidth of a signal can 
arbitrarily be changed. In addition, the linear phase remover 903 obtains an 
initial phase angle of the clock signal under measurement simultaneously with 
5 the instantaneous phase noise, and outputs the obtained initial phase angle to die 
deterministic clock skew estimator. 

A processing procedure in this timing jitter estimator 900 will be 
explained with reference to Fig. 20. In step 1001, an inputted clock signal 
under measurement is transformed by the analytic signal transformer 901 into an 

10 analytic signal whose predetermined frequency components are selectively 
passed. In step 1002, an instantaneous phase of the clock signal under 
measurement is estimated by the instantaneous phase estimator 902 using the 
analytic signal. In step 1003, a linear instantaneous phase corresponding to an 
ideal clock signal is estimated by the linear phase remover 903 from the 

1 5 instantaneous phase to obtain an initial phase angle of the clock signal under 
measurement. In step 1004, the linear instantaneous phase is removed by the 
linear phase remover 903 from the instantaneous phase to estimate an 
instantaneous phase noise A^(t). At the same time, in step 1005, timings 
(approximated zero-crossing points) closest to zero-crossing points of a real part 

20 of the analytic signal are detected by the zero-crossing detector 904 from the real 
part of the analytic signal using the previously explained zero-crossing point 
detection method. Finally, in step 1006, only the instantaneous phase noise 
data at the approximated zero-crossing points from the instantaneous phase noise 
are sampled by the zero-crossing sampler 904 to estimate a timing jitter sequence 

25 A(j>' [n] , and the process ends . 
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In the analytic signal transformer 901 used in the timing jitter estimator 
900, for example as shown in Fig. 19, only components around a fimdamental 
frequency are extracted by a band-pass filter 1101 from the clock signal under 
measurement to band-limit the clock signal under measurement, hi addition, 
5 the band-limited clock signal under measurement is inputted to a Hilbert 

transformer 1102 to Hilbert-transform this signal, and an output of the band-pass 
filter 1 101 is outputted as a real part of the analytic signal. An output of the 
Hilbert transformer 1102 is outputted as a constant part as an imaginary part o f 
the analytic signal. The band-pass filter 1101 may be either an analog filter or a 

10 digital filter, or may be implemented using a digital signal process such as FFT 
or the like. In addition, the band-pass filter 1101 may be constructed such that 
the pass bandwidth of the signal can arbitrarily be changed. 

Fig. 21 shows another configuration example of the analytic signal 
apparatus t ransformer 901 used in the timing jitter estimator 900. For example, 

15 FFT (Fast Fourier Transform) is applied to the clock signal under measurement 
by a time domain to frequency domain transformer 130 1 to transform a signal in 
time domain into a both-sided spectrum signal (for example. Fig. 11) in 
frequency domain. Negative frequency components of the transformed 
both-sided spectrum signal in frequency domain are replaced by zero by the 

20 bandwidth limiter 1302 to obtain a single-sided spectrum signal. In addition, 
regardin g, thi s singl e s ided, this single-sided spectrum signal, only components 
around the fimdamental frequency of the clock signal under measurement are 
retained and the other frequency components are replaced by zeros to band-limit 
the signal in frequency domain. Inverse FFT is applied by a frequency domain 

25 to time domain transformer 1303 to the band-limited single-sided spectrum 
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signal to transform the signal in frequency domain into an analytic signal in time 
domain. 

Fig. 22 shows further another configuration example of the analytic 
signal transformer 901 used in the timing jitter estimator 900. This analytic 
5 signal transformer 1500 comprises a buffer memory 1501 for storing therein a 
clock signal under measurement, a signal extractor rwaveform data selector) 
1502 for extracting the signal in the sequential order from the buffer memory 
1501 such that the signal being extracted is partially overlapped with the signal 
exfracted just before, a window function multiplier 1503 for multiplying each 

10 extracted partial signal by a window function, a time domain to frequency 

domain transformer 1504 for transforming each partial signal multiplied by the 
window function into a both-sided spectrum signal in frequency domain, a 
bandwidth hmiter 1505 for extracting only components around a positive 
fundamental frequency of the clock signal under measurement from the 

15 both-sided spectrum signal transformed in frequency domain, a frequency 

domain to time domain transformer 1506 for inverse-transforming an output of 
the bandwidth limiter 1505 into a signal in time domain, and an inverse window 
function multiplier 1507 for multiplying the transformed signal in time domain 
by an inverse number of the window function to obtain a band-limited analytic 

20 signal. The time domain to frequency domain transformer 1504 and the 

frequency domain to time domain transformer 1506 may be packaged using FFT 
and inverse FFT, respectively. In addition, the bandwidth limiter 1505 may be 
constructed such that the pass bandwidth of the signal can arbitrarily be changed. 
The operation in the case of fransforming the clock signal under 

25 measurement into a band-limited analytic signal using this analytic signal 
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transformer 1500 will be described with reference to Fig. 23. First, in step 
1601, the buffer memory 1501 stores therein the clock signal under 
measurement. Next, in step 1602, the signal extractor (waveform data selector) 
1 502 extracts a portion of the stored signal from the buffer memory 1501. In 
5 step 1603, the window function multiplier 1503 multiplies the extracted 
sectioned signal by the window fimction. In step 1604, the time domain to 
frequency domain transformer 1504 applies FFT to the partial signal multiplied 
by the window fiinction to transform the signal in time domain into a both-sided 
spectrum signal in frequency domain. In step 1605, the bandwidth limiter 1505 

10 replaces negative frequency components of the transformed both-sided spectrum 
signal in frequency domain with zero. Furthermore, in step 1606, the 
bandwidth limiter 1505 retains only components around the fimdamental 
frequency of the clock signal under measurement in the single-sided spectrum 
signal produced by replacing negative frequency components with zero and 

1 5 replaces the other frequency components with zero to limit the bandwidth of the 
signal in frequency domain. In step 1607, the frequency domain to time 
domain transformer 1506 applies inverse FFT to the band-limited single-sided 
spectrum signal in frequency domain to transform the signal in frequency 
domain into a signal in time domain. In step 1608, the inverse window function 

20 multiplier 1507 multiphes the inverse-transformed signal in time domain by an 
inverse number of the window fiinction used in the multiplication in the step 
1603 to obtain a band-limited analytic signal. Finally, in step 1609, a check is 
made to see if there is unprocessed data remaining in the buffer memory 1603. 
If there is unprocessed data remaining in the buffer memory, the signal extractor 

25 1502 extracts, in step 1610, the signal in the sequential order from the buffer 
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memory 1501 such that the signal being extracted is partially overlapped with 
the signal extracted just before, and thereafter the steps 1603, 1604, 1605, 1606, 
1607, 1608 and 1609 are repeated. If there is no unprocessed data, the process 
ends. The processing sequence of the step 1605 and the step 1606 may be 
5 exchanged. That is, only components around the fundamental frequency of the 
signal under measurement are retained and the other frequency components are 
replaced with zeros first to limit the bandwidth of the signal in frequency 
domain, and thereafter negative frequency components of the both-sided 
spectrum signal may be replaced with zero. 

10 In the linear phase remover 903 in the timing jitter estimator 900 shown 

in Fig. 19, for example as illustrated in the figure, an inputted instantaneous 
phase is converted into a continuous instantaneous phase by a continuous phase 
converter 91. Regarding the continuous instantaneous phase, its instantane o us 
linear phas e , i. e ., a linear instantaneous phase corresponding to a jitter-free ideal 

15 signal is estimated by a linear phase estimator 92 using, for example, a linear 
trend estimation method, i.e., by applying a linear line fitting by least squares 
method to the continuous instantaneous phase, to output an initial phase angle 
^o\t) of the clock signal under measurement Xj(t). In addition, the linear 
instantaneous phase is subtracted from the continuous instantaneous phase by a 

20 subtractor 93 to output an instantaneous phase noise A(|)(t). 

Further, Fig. 19, Fig. 21, and Fig. 22 are shown in the international 
publication WOOO/46606 (pubhshed on August 10, 2000) report. 

As indicated by dashed lines in Fig. 13, analog clock signals under 
measurement Xj(t) and Xk(t) may be digitized by AD converters 1701a and 

25 1701b, respectively to convert the analog clock signals into digital signals, and 
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the digital signals may be inputted to the respective timing jitter estimators 101a 
and lOlb. In addition, as indicated by dashed lines in Fig. 13, waveform 
clippers 1901a and 1901b may be provided to supply those input signals Xj(t) and 
Xk(t) to the AD converters 1701a and 1701b or the timing jitter estimators 101a 
5 and 101b, respectively after removing AM components in the state that phase 
modulation components, that are jitter components of the respective input 
signals, are retained in the input signals. The waveform clippers 1901a and 
1901b may be provided in the output sides of the AD converters 1701a and 
1701b, respectively. 

10 In addition, as indicated by dashed lines in Fig. 19, low frequency 

components of the instantaneous phase noise may be removed by a low 
frequency component remover 2101 from the instantaneous phase noise A(j)(t) 
outputted from the linear phase remover 903, and the instantaneous phase noise 
from which low frequency components have been removed may be supplied to 

15 the zero-crossing sampler 905. 

In the above description, the instantaneous phase noise A^{t) is sampled 
at approximated zero-crossing points to obtain a timing jitter sequence A(|)'[n]. 
However, since the linear phase remover 903 has a configuration as shown in 
Fig. 19, for example as indicated by dashed lines in Fig. 24, the sampling at 

20 approximated zero-crossing points may be inserted in series between the 
instantaneous phase estimator 902 and the continuous phase converter 91. 
Altematively, the sampling at approximated zero-crossing points may be inserted 
in series between the continuous phase converter 91 and the linear phase 
estimator 92/subtractor 93. In such a configuration, a timing jitter sequence 

25 A^[n] can also be obtained from the subtractor 93. 
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In addition, since the estimation of an instantaneous phase noise A(t)(t) 
from an instantaneous phase is performed by the configuration shown by the 
Hnear phase remover 903 in Fig. 19, the processing procedure is, as shown in 
Fig. 25, that after an instantaneous phase is obtained in the step 1002 in Fig. 20, 
5 in step 1003a, the instantaneous phase is converted into a continuous 

instantaneous phase by the continuous phase converter 91, and in step 1003b, a 
linear instantaneous phase of the continuous instantaneous phase is estimated by 
the linear phase estimator 92 from the continuous instantaneous phase. After 
that, in step 1004, the linear instantaneous phase is removed by the subtractor 93 
10 from the continuous instantaneous phase to obtain an instantaneous phase noise 

my 

Therefore, similarly to the configuration shown in Fig. 24, as indicated 
by dashed lines in Fig. 25, the approximated zero-crossing sampling may be 
applied to the instantaneous phase in step 2001 after the step 1002 to obtain a 

15 sample sequence of the instantaneous phase. And then the process may move 
to the step 1003a to transform the sample sequence into a continuous 
instantaneous phase. 

Alternatively, in step 2002, the continuous phase obtained in the step 
1003a may be sampled at the approximated zero-crossing points to obtain a 

20 sample sequence of the continuous instantaneous phase, and then the process 
may move to step 1003b to estimate an instantaneous linear phase from the 
sample sequence of the continuous instantaneous phase. In either case, in step 
1004, there is obtained a timing jitter sequence A^[n] by sampling the 
instantaneous phase noise at the approximated zero-crossing points. 

25 Regarding the deterministic clock skew estimator 102 in Fig. 14, for 
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example, as shown in Fig.26, the clock signals under measurement Xj(t) and Xk(t) 
may be inputted to zero-crossing timing detectors 81 and 82, respectively. 
Then zero-crossing timing sequences t'zero.cCn) and t'^zero.cCn) of those respective 
signals may be inputted to a subtractor 83 to obtain a time difference sequence 
5 between corresponding zero-crossing time points of the zero-crossing timing 
sequences t'zero,c(n) and t'^zero.cCn). Then a mean value of those time differences 
may be calculated by a mean value calculator 84 to obtain the mean value as a 
deterministic clock skew value r'^ew. 

The apparatus shown in Fig. 13 and Fig. 17 can be functioned by 

10 executing programs in a computer. 

According to the clock skew measurement apparatus and the clock skew 
measurement method of the present invention, a random spread (random 
components) of clock skew can be measured, and therefore an efficiency of 
clock skew test can greatly be improved. Moreover, if necessary, the apparatus 

15 can be constructed such that one of the clock skew components that is 
deterministically determined by clock distribution network paths are also 
simultaneously measured. 

In addition, according to the embodiments of the clock skew 
measurement apparatus and the clock skew measurement method of the present 

20 invention, a clock skew between clock signals each having a different frequency 
from one another can be obtained. By this, a clock signal having relatively low 
frequency can be used as a reference clock signal in a clock skew test, and 
therefore an efficiency and a practical usability of a clock skew test can greatly 
be improved. 

25 In addition, according to the embodiments of the clock skew 
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measurement apparatus and the clock skew measurement method of the present 
invention, a clock skew between clock signals under measurement can be 
estimated by simultaneously sampling in the sequential order clock waveforms 
under measurement utilizing an apparatus that can simultaneously measure two 
5 channels. By this, the number of required simultaneous samplings can be 
decreased from nC2(=N(N-1)/2) to (N-1) two-channel simultaneous 
measurements, and therefore a measurement time of clock skew can greatly be 
decreased. 

In addition, according to the clock skew measurement apparatus and the 
10 clock skew measurement method of the present invention, only the minimum 
number of pins are required when clock signals to be distributed within, for 
example, a semiconductor chip are taken out to the outside of the chip. 
Therefore, the test expense of VLSI test can be reduced. 
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WHAT IS CLAIMED IS 

1 . An apparatus for measuring a clock skew between a plurality of 
clock signals under measurement comprising: 

a timing jitter estimator to which the plurality of clock signals under 
measurement are inputted for estimating their respective timing jitter 
sequences; and 

a clock skew estimator to which the plurality of timing jitter 
sequences are inputted for calculating a timing difference sequence between 
those timing jitter sequences to output a clock skew sequence. 

2. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 1 further including: 

a second clock skew estimator to which a plurality of the clock skew 
sequences are inputted for obtaining a difference between the plurality of 
clock skew sequences. 

3 The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 2 further including: 

a frequency multiplier to which the timing jitter sequence is inputted 
for multiplying a frequency of the timing jitter sequence to output the 
frequency-multiplied timing jitter sequence to said clock skew estimator. 

4 The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 1 further including: 

a frequency multiplier to which the timing jitter sequence is inputted 
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for multiplying a frequency of the timing jitter sequence to output the 
frequency-multiplied timing jitter sequence to said clock skew estimator. 

5. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 1 further including a 
deterministic clock skew estimator for estimating a timing error between 
ideal clock edges of the plurality of clock signals under measurement to 
output a deterministic component of clock skew to said clock skew 
estimator, wherein said clock skew estimator is an estimator for adding the 
deterministic component of clock skew to the timing difference sequence to 
output the summed value as the clock skew sequence. 

6. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to any one of claims 1-5 further 
including: 

a clock skew detector to which the clock skew sequence is inputted 
for obtaining clock skew values of the clock signals under measurement 
from the clock skew sequence. 

7. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 6, wherein said clock 
skew detector comprises one or a plurality of a peak-to-peak detector for 
obtaining a difference between the maximum value and the minimum value 
of the clock skew sequence, an RMS detector for obtaining a 
root-mean-square value of the clock skew sequence, and a histogram 
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estimator for obtaining a histogram of the clock skew sequence. 

8. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to any one of claims 1-5, 
wherein said timing jitter estimator comprises; 

an analj^ic signal transformer for transforming a clock signal under 
measurement into a complex analytic signal; 

an instantaneous phase estimator for obtaining an instantaneous 
phase of the analytic signal; 

a continuous phase converter for converting the instantaneous phase 
into a continuous instantaneous phase; 

a linear phase estimator for estimating, from the continuous 
instantaneous phase, its linear instantaneous phase; 

a subtractor for removing the linear instantaneous phase from the 
continuous instantaneous phase to obtain an instantaneous phase noise; and 

a zero-crossing sampler for sampling its input at timings close to 
zero-crossing timings of a real part of the analytic signal to output the 
sampled signal, said zero-crossing sampler being inserted in series to any 
one of connection points between said instantaneous phase estimator and 
said continuous phase converter, between said continuous phase converter 
and said linear phase estimator/subtractor, and at an output side of said 
subtractor; 

wherein a timing jitter sequence of the clock signal under 
measurement is outputted as an output o f said jitter sequence estimation part. 
said timing jitter estimator . 
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9. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 8, wherein said 
deterministic clock skew estimator of the plurality of clock signals under 
measurement is an estimator that obtains a difference between initial phase 
angles of the linear instantaneous phases to obtain a deterministic 
component of clock skew. 

10. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 8, wherein said analytic 
signal transformer can change a pass bandwidth of the clock signal under 
measurement. 

1 1 . The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to claim 8, wherein said timing 
jitter estimator includes a low frequency component remover to which the 
instantaneous phase noise is inputted for removing low frequency 
components of the instantaneous phase noise to output the instantaneous 
phase noise from which the low frequency components have been removed. 

12. The apparatus for measuring a clock skew between a plurality of 
clock signals under measurement according to any one of claims 1-5 further 
including; 

a waveform clipper to which the clock signal under measurement is 
inputted for removing amplitude modulation components of the clock signal 
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under measurement in the state that phase modulation components are 
retained in the clock signal under measurement to output the clock signal 
under measurement from which the amplitude modulation components have 
been removed. 

13 . A method of measuring a clock skew between a plurality of clock 
signals imder measurement comprising: 

a step of estimating timing jitter sequences of the respective clock 
signals under measurement; and 

a step of calculating a timing difference sequence between the 
plurality of timing jitter sequences to estimate a clock skew sequence. 

14. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 13 further including: 

a step of obtaining a difference between the plurality of clock skew 
sequences to estimate a clock skew sequence. 

15. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 14 further including: 

a step of assign each timing jitter of the timing jitter sequence M 
times to estimate a timing jitter sequence that is created by multiplying a 
frequency of the corresponding clock signal under measurement by (M+1). 

16. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 13 further including: 

a step of assign each timing jitter of the timing jitter sequence M 
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times to estimate a timing jitter sequence that is created by multiplying a 
frequency of the corresponding clock signal under measurement by (M+1). 

17. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 13 further including a step of 
estimating a timing error between ideal clock edges of the plurality of clock 
signals under measurement to estimate a deterministic component of clock 
skew, wherein said step of estimating a clock skew sequence is a step of 
adding the deterministic component of clock skew to the timing difference 
sequence to obtain the summed value as the clock skew sequence. 

18. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to any one of claims 13-17 further 
including: 

a step of obtaining clock skew values of the clock signals under 
measurement from the clock skew sequence. 

19. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 18, wherein said step of 
obtaining a clock skew comprises one or a plurality of a step of obtaining a 
difference between the maximum value and the minimum value of the clock 
skew sequence to calculate a peak-to-peak value, a step of obtaining a 
root-mean-square value of the clock skew sequence to calculate an RMS 
value, and a step of obtaining a histogram data of the clock skew sequence. 
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20. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to any one of claims 13-17, wherein 
said step of estimating a timing jitter sequence comprises; 

a step of transforming a clock signal under measurement into a 
complex analytic signal; 

a step of obtaining an instantaneous phase of the clock signal under 
measurement from the analytic signal; 

a step of converting the instantaneous phase into a continuous 
instantaneous phase; 

a step of estimating, from the continuous instantaneous phase, its 
linear instantaneous phase; 

a step of removing the linear instantaneous phase from the 
continuous instantaneous phase to obtain an instantaneous phase noise; and 

a step of sampling any one of the instantaneous phase, the continuous 
instantaneous phase, and the phase noise waveform at timings close to 
zero-crossing timings of a real part of the analytic signal, 

wherein a timing jitter sequence of the clock signal under 
measurement is eventually obtained. 

21 . The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 20, wherein said step of 
estimating a deterministic component of clock skew between the clock 
signals under measurement is a step of obtaining a difference between initial 
phase angles of linear instantaneous phases of the plurality of clock signals 
under measurement to obtain a deterministic component of clock skew. 
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22. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 21, wherein said step of 
estimating a deterministic component of clock skew between the clock 
signals under measurement comprises: 

a step of obtaining an offset signal in which either a correlation 
between timing jitter sequences of the plurality of clock signals under 
measurement or a correlation between instantaneous phase noises of the 
plurality of clock signals under measurement shows the maximum value to 
obtain an offset value of clock edge; and 

a step of obtaining a sum of the offset value and the difference 
between the initial phase angles to obtain the deterministic component of 
clock skew. 

23. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 20, wherein said step of 
estimating a deterministic component of clock skew between the clock 
signals under measurement is a step of obtaining a mean value of differences 
of zero-crossing timings between the plurality of signals under measurement 
to obtain a deterministic component of clock skew. 

24. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to claim 20, wherein said step of 
estimating a timing jitter includes a step of removing low frequency 
components of the instantaneous phase noise. 
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25. The method of measuring a clock skew between a plurality of clock 
signals under measurement according to any one of claims 13-17 further 
including; 

a step of performing a waveform clipping in the state that phase 
modulation components of the clock signal under measurement are retained 
to remove amplitude modulation components of the clock signal under 
measurement, and moving to the step of estimating a timing jitter sequence. 
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ABSTRACT 



Timing jitter sequences A(|)^[n] and A^\n] of respective clock signals 
under measurement Xj(t) and Xk(t) are estimated, and a timing difference 
sequence between those timing jitter sequences is calculated. In addition, 
initial phase angles and of linear instantaneous phases of the Xj(t) and 
Xk(t) are estimated, respectively. A sum of a difference between those 
initial angles and the timing difference sequence is calculated to obtain a 
clock skew sequence between the Xj(t) and Xk(t). 



